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ABSTRACT 

The properties of the first generation of black-hole seeds trace and distinguish different models of 
formation of cosmic structure in the high-redshift universe. The observational challenge lies in identifying 
black holes in the mass range ~ 100—1000 M Q at redshift z ~ 10. The typical frequencies of gravitational 
waves produced by the coalescence of the first generation of light seed black-hole binaries fall in the 
gap between the spectral ranges of low-frequency space-borne detectors (e.g., LISA) and high-frequency 
ground-based detectors (e.g., LIGO, Virgo and GEO 600). As such, these sources are targets for proposed 
third-generation ground-based instruments, such as the Einstein Telescope which is currently in design 
study. Using galaxy merger trees and four different models of black hole accretion — which are meant 
to illustrate the potential of this new type of source rather than to yield precise event-rate predictions 
— we find that such detectors could observe a few to a few tens of seed black-hole merger events in three 
years and provide, possibly unique, information on the evolution of structure in the corresponding era. 
We show further that a network of detectors may be able to measure the luminosity distance to sources 
to a precision of ~ 30%, allowing us to be confident of the high-redshift nature of the sources. 

Subject headings: black hole physics — gravitational waves 



to ~ 10 3 M© with a signal-to-noise ratio (SNR) ;> 10, but 
will not probe lower masses. To probe the 10 2 — 10 3 M Q 
range, a new GW telescope with optimal sensitivity in the 
frequency window 0.1 — 10 Hz is needed. This could be 
achieved either by a second generation space-based detec- 
tor, such as the Big Bang Observer (Phinney et al. 2003), 
ALIA (Bender et al. 2005) or DECIGO (Kawamura et 
al. 2006), or by the third generation of ground-based laser 
interferometers, for which the target is ~ ten-fold strain 
sensitivity increase over Advanced LIGO and a low fre- 
quency cut-off at ~ 1 Hz (the ET target sensitivity is 
compared to other instruments in Fig. 1). In this Letter, 
we focus on third-generation ground-based instruments as 
these are presently undergoing conceptual design studies, 
and for reference we use the specific example of the Ein- 
stein gravitational- wave Telescope (ET) (see e.g., Freise et 
al 2009). We demonstrate that instruments such as the ET 
can detect seed black-hole binaries, albeit at likely rates 
of no more than a handful a year, and discuss whether 
such observations can uniquely identify these events as 
produced by light remnants of Pop III stars. 

2. MODELS FOR POP III SEED GROWTH 

The astrophysical scenario that we consider in this Let- 
ter assumes that seed black holes of mass ~ 100 M Q are 
produced by the first generation of supernovae in the very 
high-redshift Universe at z as 20. These black holes are 
efficient at accreting mass and hierarchically merge fol- 
lowing mergers between their host halos. We can trace 
the merger hierarchy by means of Monte-Carlo merger-tree 
realizations based on the extended Press-Schechter formal- 
ism (Press & Schcchtcr 1974), assuming a standard ACDM 

1 Center for Gravitational Wave Physics, The Pennsylvania State University, University Park, PA16802, USA 
institute of Astronomy, University of Cambridge, Cambridge, CB3 0HA, UK 
3 Department of Physics and Astronomy, Northwestern University, Evanston, IL, 60208 USA 
4 School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK 

1 



1. INTRODUCTION 

An ability to probe the nature of the first massive black- 
hole (MBH) seeds at medium-to-high redshift is fundamen- 
tal to understanding the hierarchical assembly scenario, 
to discriminate between different models of structure for- 
mation in the high-redshift universe and to explore the 
link between black holes residing at the centre of galax- 
ies and their hosts (Sesana, Volonteri & Haardt 2007). If 
MBH seeds are massive (i.e., ~ 10 5 M©), the future space- 
borne gravitational wave (GW) detector LISA (Bender et 
al. 1998) will probe the first epoch of mergers between 
these seeds at high redshift. However, if MBH seeds arc 
light (i.e., ~ 100 M Q ), the GWs from these mergers will 
fall between the sensitive frequency band of LISA and 
of currently operating and planned Advanced versions of 
ground-based instruments — LIGO, Virgo and GEO 600 
(see Sigg et al. 2008, Acernese et al. 2008 and Grote et 
al. 2008 for recent status reports). In this Letter we show 
that third-generation laser interferometers may be able to 
fill this gap by directly probing the first mergers between 
light MBH seeds and will thus provide complementary in- 
formation to other instruments. 

If seed black holes (BHs) are the remnants of Pop III 
stars with mass ~ 100 M© (Madau & Rees 2001; see 
Sesana, Volonteri & Haardt 2007 for a short review of 
this and alternative scenarios) we expect dozens of MBH 
binary (MBHB) coalescences per year in the mass range 
- 10 2 - 10 6 M (Sesana et al. 2004) (NB we will use 
MBHB liberally to refer to any binary formed between 
black holes in the centres of merging dark-matter halos). 
Most of the MBHB events occur at redshift z ;> 3, with the 
consequence that LISA will be able to detect MBHBs down 
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cosmology with the 1-year WMAP cosmological parame- 
ters (Spergel et al. 2003); technical details are given in 
Volonteri, Haardt & Madau (2003; VHM) and Volonteri, 
Salvaterra & Haardt (2006; VSH). We consider four vari- 
ants of this model, which have the same merger history for 
the dark-matters halos, but different prescriptions for the 
mass-distribution and accretion efficiency of the seeds: (i) 
VHM, ems (VHM with equal-mass seeds) is based on equal 
150 Mq seeds accreting at the Eddington limit during each 
merger episode; (ii) VHM.smd (VHM with seed-mass dis- 
tribution) differs from VHM, ems only in the seed-mass dis- 
tribution, which is log-normal in the range 10 — 600 M Q ; 
(iii) Shank has black-hole seeds with flat mass distribution 
in the range 150 — 600 M Q and redshift- dependent accre- 
tion efficiency, following Shankar et al. 2004; (iv) Hopk 
again assumes black-hole seeds with flat mass-distribution 
in the range 150 — 600 M and accretion efficiency that is 
luminosity-dependent according to the prescription given 
by Hopkins et al. 2005. 

Integrated over cosmic history, all these models repro- 
duce the X-ray and optical quasar luminosity function at 
z < 3, the observed faint X-ray counts of AGNs (see 
VSH), and, integrated over all MBH masses, predict about 
50 — 70 black-hole coalescences per year in the Universe. 
This range is statistical and does not include uncertainties 
in assumptions about, e.g., cosmology, that go into the 
merger trees, which could change the number of predicted 
events by a factor of a few in either direction. The accre- 
tion models considered here have been tuned to reproduce 
observations at low redshift, z <J 3, rather than accretion 
onto seed black holes at z ;> 5. Recent work has indi- 
cated that accretion onto IOOMq black holes is inefficient 
(Alvarez et al. 2008, Milosavljevic et al. 2008) and gen- 
erally sub-Eddington. This adds further uncertainties to 
the light seed scenario, but we emphasize that our choice 
of models is guided by the goal of illustrating the science 
potential of this new class of observations rather than an 
attempt to provide solid predictions for event rates. 
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Fig. 1. — The thin curves show the characteristic amplitudes fh(f) 
of the frequency-domain gravitational waveforms IMR (solid red) 
and EOBNR (dashed blue) for an equal-mass, optimally-oriented 
source with a redshifted total mass of 5OOM0 at a luminosity dis- 
tance of 33 Gpc. The thick curves show the characteristic noise 

amplitude spectral densities \J fS n {f) for LISA (dotted magenta), 
Enhanced LIGO (solid cyan), Advanced LIGO (dashed green) and 
ET (dash-dotted black). The ratio of the characteristic amplitudes 
of the waveform to the noise yields the SNR accumulated in a band- 
width equal to the frequency. 



3. WAVEFORM AND DETECTOR MODELS 

To compute the sensitivity of the ET to GWs generated 
during the coalescence of seed MBH binaries, we model 
the gravitational wave emission with the phcnomcnological 
waveform family (IMR) introduced by Ajith et al. (2008). 
This describes radiation from a non-spinning black-hole bi- 
nary and includes in a self-consistent manner the inspiral, 
merger and ring-down phases; these phases are marked in 
Fig. 1, which shows a typical frequency-domain gravita- 
tional waveform h(f) = A{f)e l ^^\ Fig. 1 also indicates 
that a significant fraction of the signal is contributed by 
the merger and ringdown, so an inspiral-only waveform 
would be inadequate. 

Exact expressions for the amplitude A(f) and the phase 
ip(f) of the IMR waveforms are provided by Eqs. (4.17), 
(4.18) and (4.19) and Table I of Ajith et al. (2008). The 
strain at the detector depends on the total redshifted mass 
M z = (1 + z)M, the symmetric mass ratio n, a fiducial 
time of arrival parameter to, and six extrinsic parameters 
- two sky-position angles, the orbital phase at time to, 
the wave polarization angle tjj, the source inclination angle 
with respect to the line of sight t, and the luminosity dis- 
tance to the source D^. The optimal signal-to-noise ratio 
(SNR) at the instrument is SNR 2 = 4 j °° df\h(f)\ 2 /S„(f), 
where S n (f) is the one-sided noise power spectral density 
of the interferometer, as shown in Fig. 1. For a network 
of detectors, the total network SNR is obtained by adding 
the SNRs of the individual instruments in quadrature. We 
report SNRs averaged over the sky location of the source 
and its orbital-plane orientation. 

As a check, we also computed the SNR for some 
events using the effective-one-body-numerical-relativity 
(EOBNR) waveform family introduced by Buonanno et 
al. (2007). The SNRs predicted by the EOBNR wave- 
forms are typically somewhat higher (by up to ~ 25%) for 
equal-mass events, and are comparable for all events ex- 
cept those with very asymmetric mass ratios, r\ <C 1, where 
neither waveform family has been shown to be valid. 

The currently favoured design for the Einstein Telescope 
calls for a 10km triangular facility containing three 60° 
detectors (Freise et al. 2009); we refer to this design as 
a "single ET". We use the noise power spectral density 
defined in Hild et al. (2008) for a single right-angle 10km 
detector. The angle-averaged sensitivities achieved by two 
right-angle detectors and three 60° detectors (single ET) 
are factors of y/2 and 3/2 higher, respectively, than for one 
right-angle 10km instrument. 

4. DETECTION-RATE ESTIMATES 

For each of the four scenarios described in Section 2, we 
have generated 1000 independent realisations of the galaxy 
merger history and computed the SNR of all coalescences 
that take place over a time span of 3 years, representative 
of a typical data-taking period for a third-generation in- 
strument. Figure 2 shows the number of events that would 
be seen by the third generation network in three years as a 
function of the SNR in a single right-angle 10km detector 
within the network, for each of the four models. Assuming 
that a network SNR of 8 would be required for detection, 
the SNR in one 10km detector would be 5.3 for a single 
ET or 4.4 (3.9, 3.1) if one (two, three) additional 10km de- 
tector (s) were added. Figure 2 indicates that in all models 
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except Shank, we could expect to detect a few to a few 
tens of events over three years. The Shank model predicts 
fewer than one event in a single ET. 




2 4 6 8 10 12 14 
SNR threshold 

Fig. 2. — Number of events detected in three years as a function of 
threshold SNR required in a single right-angle 10km detector. For 
a network SNR of 8, this threshold would be 5.3 for a single ET, or 
4.4 (3.9, 3.1) with one (two, three) additional 10km detector(s). 

The ET will be sited underground, which should make it 
possible to suppress seismic gravity-gradient noise (Hughes 
& Thorne 1998) and achieve sensitivity at frequencies as 
low as ~ 1 Hz (cf. Advanced LIGO, Fig. 1); however, it 
is not currently clear whether it will be possible to mit- 
igate this noise source below ~ 3 Hz. In addition, the 
ET may suffer a confusion background in the l-5Hz range 
due to cosmological compact-binary systems (Regimbau & 
Hughes 2009). The impact of such a background on Pop 
III event rates needs to be properly quantified, but we have 
checked this crudely. The preceding results assumed a low- 
frequency cut-off of 1Hz (which is also the cut-off we use 
for parameter estimation), but we found that only ~ 25% 
of the events were lost with a higher cut-off at 5Hz. We 
conclude that design changes or a confusion background 
should not affect our qualitative conclusions. 

Figure 3 indicates the distribution of masses and red- 
shifts for the events that would be detected under each 
scenario. We see that most events have intrinsic masses of 
a few hundred Mq (down to a few tens in the VHM,smd 
model), although a few events with M ^> 1000 M Q might 
also be observed. Many of the events will be at medium- 
to-high redshift, 5 <J z <, 10, except in the VHM,smd 
model which predicts a distribution peaked around z ~ 4 
with a long tail extending to z > 12. Most of the ob- 
servable events are related to the hierarchical assembly of 
small-to-medium size parent halos (masses in the range 
10 11 — 10 13 M Q at z = 0). In such halos, seeds are less 
likely to experience coalescences accompanied by major ac- 
cretion episodes at high redshift, and they arc more likely 
to have M <J 1000 M Q at z = 10, leaving them in the band 
accessible to third generation detectors. 

The number of events and their redshift distribution 
may provide constraints on both the mass function and 
accretion history of seed black holes in the early Universe. 
Considering our example scenarios, in the Shank model 
seeds are born quite massive (150 — 600 Mq) and they ac- 
crete at the Eddington limit at high redshift [see equation 
(2) in VSH]; in this model, most of the seeds grow above 



1000 M Q by z = 10, making detection with ground-based 
detectors difficult. In the Hopk model, accretion is on av- 
erage less efficient and many of the seeds are still in the 
few- hundred solar-mass range at z < 10, making detec- 
tion easier. The initial seed mass distribution may also 
leave a detectable signature on the observed events. The 
VHM,ems and VHM,smd models are characterised by the 
same Eddington-limited accretion prescription; nonethe- 
less, their mass and redshift distributions are significantly 
different. In the VHM,ems model most of the 150 M Q 
seeds grow at high z, radiating outside the l-1000Hz band; 
only the 'tail' of <; 10 3 M black holes left behind at z < 10 
is detectable. In the VHM,smd model, there are many 
< 100 Mq seeds which would be observable at very high 
redshift (z > 10). Many of them will grow inefficiently, 
and still have masses ~ 100 M Q at low redshift, making 
them perfect targets for ET; however, it may be impossible 
to discriminate between coalescences involving these low- 
z Pop III remnants and those involving intermediate-mass 
BHs (IMBHs) formed at low redshift. IMBHs might form 
via runaway collisions of massive stars in globular clusters 
(Portegies-Zwart ct al. 1999; but see Glebbeek et al. 2009), 
and scenarios have been proposed in which IMBH binaries 
also form (Fregeau et al. 2006). The IMBH masses consid- 
ered there, ~ 10 3 Mq, are somewhat higher than seed black 
hole masses, ~fewxl00M Q . Therefore, if we observe an 
IMBH merger at z J> 5 with masses ~ 1OOM , we can be 
fairly confident that the constituent black holes are Pop 
III seeds. Further work will be needed to identify what 
observational signature (s) provide the best discriminating 
power between these two formation channels. 
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Fig. 3. — Distribution of source-frame masses (top) and redshifts 
(bottom) for detected (SNR> 5) events, normalised such that the 
integral over the distribution is unity. 
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This current work indicates that, while probing seeds 
mostly at z < 10, third-generation detectors may provide 
useful insights into the properties and accretion history of 
seed black holes at higher redshift. LISA will be able to 
probe deeper in redshift, since ground-based detectors are 
limited at low frequency by a gravity-gradient 'noise wall' 
at 1Hz. However, LISA will not be able to observe seed 
black-hole mergers because it is not sufficiently sensitive in 
the relevant frequency range (see Fig. 1). LISA will only 
provide indirect constraints on seed black-hole populations 
through observations of subsequent mergers of black holes 
in the mass-range 10 4 -10 6 M@. 

5. PARAMETER-ESTIMATION RESULTS 

In order to identify an event as a seed black-hole merger, 
we must be able to determine the redshift and mass of 
the system. As mentioned earlier, events at low red- 
shift z <J 3 might not involve primordial black holes. A 
gravitational-wave observation determines the redshifted 
total mass (1 + z) M and luminosity distance Dl{z), but 
not the redshift z independently. It is unlikely that electro- 
magnetic counterparts to ET events will be observed, but 
if we assume a concordance cosmology inferred from other 
observations, we can determine the redshift from the lumi- 
nosity distance. The fractional redshift error is then com- 
parable to the fractional luminosity-distance error, plus an 
0(10%) error from uncertainties in the cosmology. 

We used the Fisher information matrix to evaluate the 
parameter-estimation errors. We carried out a Monte 
Carlo simulation over the extrinsic parameters for fixed 
choices of the intrinsic parameters, M z and rj. We note 
that this approach may overstate the precision of parame- 
ter estimation at low SNRs (Vallisneri 2008); a more rigor- 
ous study of parameter-estimation accuracy would require 
the computation of the full posterior probability density 
function and is beyond the scope of the present paper. 

The intrinsic parameters, M z and 77, are determined by 
the waveform phase evolution and so we expect to measure 
them well. The extrinsic parameters, by contrast, affect 
only the relative amplitudes of the signal. A single ET 
makes four independent measurements (two quadratures 
in two detectors) and so at least one additional interfer- 
ometer is needed to determine the six extrinsic parameters. 
We computed parameter-estimation errors for several net- 
work configurations consisting of an ET at the location of 
Virgo and additional interferometers at Hanford and Liv- 
ingston. At a fixed network SNR of 8, we found that for 
all networks M z and 77 would be determined to a fractional 
accuracy better than 1% for all but the most massive sys- 
tems (M z ~ 1000A/ Q ). The distance, and hence redshift, 
will be determined less accurately. One ET plus a second 
10km interferometer will be enough to determine the dis- 
tance to better than ~ 40% in 68% of cases. The addition 
of a third 10km detector will improve this distance preci- 
sion to ~ 30%. Upgrading the 10km detectors to ETs leads 
to further modest improvements, but would also increase 
the SNR accumulated from a source at a fixed distance. 
We also explored siting the second detector in Perth, Aus- 
tralia, instead of Hanford, but found little net difference in 
parameter-estimation accuracy. The error in the source- 
frame mass is dominated by the error in z rather than that 
in M z , and so is also ~ 40%. Therefore, we should be able 
to say confidently that a source at z ~ 5 is at high red- 



shift and has M ~ 100 M , and hence is most likely a seed 
black-hole merger. 

6. DISCUSSION 

Information about seed black holes is extremely difficult 
to obtain using present or future electromagnetic observa- 
tions. The only direct means to study seed black holes 
is via gravitational-wave observations of mergers. If seed 
black holes are remnants of Pop III stars at high redshift, 
then LISA may not be able to probe the early stages of 
their evolution. In this paper we have analyzed the ability 
of third-generation ground-based interferometers, such as 
the proposed Einstein Telescope, to detect the coalescences 
of 100 — 1000 M Q seed black-hole binaries for a range of 
seed properties and accretion histories. 

We have found that third-generation detectors will be 
able to detect ~ 1-30 events over a 3-year observation, 
depending on the selected model and on the assumed tele- 
scope configuration. The distribution of detected masses, 
ranging from a fewxl0M Q to a fewx 1000 M , is comple- 
mentary to the range probed by LISA, making the detec- 
tion of low-mass seeds possible. The noise wall at 1 Hz will 
preclude the detection of sources at very high redshift, al- 
though in the case of a seed-mass spectrum that extends 
down to 10 M Q (i.e., the VHM,smd model), some detec- 
tions may be possible at z > 10. The mass and redshift 
distribution of detected events may be useful in recon- 
structing the accretion history of the first seeds. We have 
also shown that a detector network with at least two sites 
will be able to determine the (redshifted) masses and lu- 
minosity distances of the majority of events to accuracies 
of <; 1% and <; 40% respectively. This should be sufficient 
for us to say with confidence that the merger is occurring 
at high redshift between low-mass seeds. 

Previous work (Wyithe & Loeb 2004) indicated that 
Adv. LIGO might detect pop III seeds. That work used a 
semi-analytic model for structure growth and our current 
work contradicts this, suggesting that their model signifi- 
cantly over-predicted the number of mergers occurring at 
low redshift. There are various sources of uncertainty in 
our analysis. Our waveform model ignored the effects of 
black-hole spin and higher gravitational-wave harmonics, 
which tend to enhance the SNR and increase the event 
rate. These effects also help to break degeneracies between 
parameters, improving the validity of the Fisher Matrix 
approach to computing parameter accuracies used in this 
paper. However, the SNR threshold of 8 may be optimistic 
when source confusion and realistic instrumental noise are 
taken into account. The merger-tree models have various 
uncertainties, such as the choice of cosmological parame- 
ters, e.g., which could change the rates by a factor of a 
few in either direction. More work is also needed to under- 
stand how to distinguish IMBH binaries formed in globular 
clusters from seed black holes. Despite these uncertainties, 
this Letter is a proof of concept which clearly demonstrates 
that third-generation ground-based instruments have the 
capability to detect seed black holes, allowing us to con- 
firm (or discard) the hypothesis that MBH seeds are light 
remnants of Pop III stars, and suggests that this could be 
one of the science drivers for these instruments. 
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